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Calorons and BPS monopoles with non-trivial holonomy in the 
confinement phase of SU(2) gluodynamics* 
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With the help of the cooling method applied to SU(2) lattice gauge theory at non-zero T < Tc we present 
numerical evidence for the existence of superpositions of Kraan-van Baal caloron (or BPS monopole pair) solutions 
with non-trivial holonomy, which might constitute an important contribution to the semi-classical approximation 
of the partition function. 
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A large number of applications to hadronic 
phenomenology has proven the instanton liquid 
model to be a powerful non-perturbative calcula- 
tional framework of QCD. Whereas fundamental 
non-perturbative features like spontaneous chiral 
symmetry breaking and the large ry'-mass can be 
understood in terms of instantons, quark con- 
finement seems to be related to other excita- 
tions like Abelian monopoles or center vortices. 
That instantons constitute important contribu- 
tions to the Euclidean QCD vacuum transition 
amplitude has found support from lattice inves- 
tigations. What possibilities exist to extend the 
present instanton model? On one hand, correla- 
tions among instantons have to be taken into ac- 
count. On the other, the question arises, whether 
other extrema of the Euclidean action also have to 
be taken into consideration in the semi-classical 
ansatz. What role do play extended monopoles? 

For T 7^ the semi-classical approach [Q so 
far was based on superpositions of time-periodic 
instantons, the so-called calorons [||. These HS 
calorons have trivial holonomy, i.e. the Polyakov 
loop at spatial infinity takes values in the cen- 
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ter of the gauge group. Recently, Kraan and 
van Baal (KvB) have constructed a more gen- 
eral class of periodic and stable solutions of the 
Euclidean Yang-Mills field equations exhibiting 
non- trivial holonomy ||^. The HS calorons turn 
out to be a limiting case of KvB solutions. The 
main feature of the SU{Nc) KvB solutions, hav- 
ing topological charge Qt = ±1 , is that they 
can dissolve into Nc separate static lumps of non- 
integer topological charge, each of them consti- 
tuting a BPS monopole. Taken as a background 
field for the Dirac fermion operator, they lead 
to zero-modes concentrated only at one of these 
BPS monopoles depending on the boundary 
conditions used for the fermion field. These fea- 
tures can be used to identify KvB solutions in 
smoothed lattice fields. 

In this contribution we are discussing some new 
results (see Q ) of continuing lattice investigations 
presented also at previous LATTICE conferences. 
For pure Yang-Mills theory at non-zero tempera- 
ture we identify low-action semi-classical configu- 
rations which are the background fields for equi- 
librium gauge fields in the path integral represen- 
tation of the partition function. We start from 
Monte Carlo (MC) generated lattice gauge fields. 
By standard relaxation ('cooling') we minimize 
the action down to metastable plateaus shown 
to be approximate solutions of the lattice field 
equations. We restrict ourselves to the SU{2) 
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case and use the Wilson action. The lattice size 
throughout this paper is 16^ x 4, i.e. the de- 
confinement transition occurs at /3 = 4/(/^ = 
Pc — 2.29 . We shall concentrate on the confine- 
ment phase. In contrast to our previous investiga- 
tions we employ standard periodic boundary con- 
ditions (p.b.c.) both for the MC simulation and 
for cooling. For the action as a function of cool- 
ing iterations we observe shoulders or plateaus 
close to multiple values of the one-instanton ac- 
tion value S'inst = 27r^/3. We stop cooling on 
plateaus when the second derivative of the ac- 
tion changes its sign from positive to negative 
values. At these stages of cooling the gauge fields 
are investigated with respect to the action den- 
sity, the naively discretized topological density 

the Polyakov loop variable L{x) — ^tr Ylu^i ^K,t,4 
and the eigenvalues A and eigenmode densities 
tpijj'^ [x) of the non-Hermitean Wilson-Dirac oper- 
ator D[U], with Y.yD[U]^.y ipiy) = A We 
have applied the Arnoldi method to both cases 
of time-antiperiodic and time-periodic boundary 
conditions of the fermion fields. The statistics for 
our investigations were 0(200) per /3-value. 

First we investigate the plateau configurations 
at S* ~ "S'inst and below. There we have seen three 
types of non-trivial objects. First of all, with the 
highest frequency of 60 to 70%, we have found 
approximate solutions which show all indications 
for (anti)selfdual KvB solutions consisting of one 
or two separate lumps of topological charge. Both 
cases are characterized by opposite-sign peaks of 
the Polyakov loop and non-trivial ('asymptotic') 
holonomy values. The Wilson-Dirac operator ex- 
hibits always one distinct real ('zero'-) mode. For 
time-antiperiodic b.c. the related eigenfunction 
is localized just at the peak of the Polyakov loop 
which corresponds to so-called Taubes winding 
(see 1^). Repeating the eigenmode investigation 
with periodic b.c. provides eigenfunctions local- 
ized at the complementary Polyakov loop peak in 
exact correspondence with the analytic results . 
We have denoted these configurations as calorons 
{CAL), if they were consisting of one topological 
lump and were non-static in the time-direction, 
whereas we called them dyon-dyon pairs {DD) 



if a dissociation into two lumps of (non-integer) 
topological charges became visible and their fields 
were static in time. A typical DD configuration 



X 1 O 




Figure 1. 2D cuts for a DD pair after cooling: 
(a) topological charge density, (b) Polyakov loop, 
(c) 'zero'-mode density for time-antiperiodic b.c. 

has been drawn in Fig. |l|. Fits of DD and CAL 
configurations with the analytic KvB solutions 
were done and have shown to work well. There- 
fore, we conclude that KvB solutions at lowest in- 
stanton action plateaus are found dominating for 
T < TcJ3 The second class of very stable config- 
urations (around 10 to 20%) were dyon-antidyon 
pairs {DD) consisting of lumps of opposite half- 
integer topological charge and of equal-sign peaks 
of the Polyakov loop variable. The fermionic 
spectrum showed pairs of almost real eigenvalues 
complex conjugate to each other. These objects 
are not described by KvB solutions and look like 
peculiar superpositions of BPS monopoles of op- 
posite topological charge. Their role is still un- 
clear. Finally, in quite rare cases (less than 10%) 
we have seen pure magnetic objects like Dirac 

^The situation is completely different for T > Tc- Self- 
dual semi-classical fields become strongly suppressed un- 
der cooling if p.b.c. are applied. 
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sheets or pure magnetic monopoles. In contrast, 
such objects seem to dominate the semi-classical 
structure at T > @. 

Higher action plateaus S ~ (2 • ■ • SjS'inst have 
been investigated for T < Tc with the same 
methods and observables. One can easily iden- 
tify superpositions of CAL's, D's and D's. The 
Polyakov loop variable as well as the fermionic 
modes with both kinds of boundary conditions 
allowed us to prove that these configurations al- 
ways exhibit a pair structure (of Polyakov line 
peaks) supporting the view of superpositions of 
KvB solutions with an admixture of DD-pairs. 
Of course, individual pairs are hardly identified. 

What is the typical holonomy ascribed 
to these configurations? We have deter- 
mined the distribution P(Loo) of the Polyakov 
loop averaged over low action regions M = 

{ ^ I W7 ^(^' — 10"^ I in the form Loo = 
Esgm L{x). The result is seen in Fig. |^. The 
histogram peaks around zero indicating the dy- 
namical dominance of non-trivial holonomy. The 
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The latter provide a completely different picture. 
Therefore, the conditional distributions -P[L|'^] 
are a hint in favor of KvB solutions, also for en- 
sembles with several lumps of topological charge. 




(b) 

Figure 3. Conditional distributions L'[i|<;] (a) at 
plateaus S ~ 4S'inst, (b) for randomly distributed 
KvB solutions. 

We would like to conclude that an improvement 
of the semi-classical caloron approach to the path 
integral at < T < T^, should take into account 
superpositions of solutions with non-trivial holon- 
omy. 



Figure 2. Histogram P(Loo) for S ~ 

4S'inst plateaus compared with the Haar measure 
(dashed fine). 

Polyakov line behaviour in different intervals of 
local action can be used in order to see that also 
complicated superpositions of D's and D's are 
similar to KvB-solutions. For low-ioo configu- 
rations we have found a local correlation between 
the Polyakov loop and the action density val- 
ues <;{x) = s{x, t) in the form of conditional 
distributions P[L\(,\. In Fig. |^ we compare the re- 
sult seen on four-instanton plateaus with that of 
lattice-discretized KvB solutions with randomly 
distributed parameters. 

We compared also with randomly distributed 
'old-fashioned' HS calorons of trivial holonomy. 
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